Phase separation dynamics in colloid-polymer mixtures: the effect of interaction range 
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Colloid-polymer mixtures may undergo either fluid-fluid phase separation or gelation. This depends on the 
depth of the quench (polymer concentration) and polymer-colloid size ratio. We present a real-space study 
of dynamics in phase separating colloid-polymer mixtures with medium- to long-range attractions (polymer- 
colloid size ratio #r = 0.45 — 0.89), with the aim of understanding the mechanism of gelation as the range of 
the attraction is changed. In contrast to previous studies of short-range attractive systems, where gelation occurs 
shortly after crossing the equilibrium phase boundary, we find a substantial region of fluid-fluid phase separation. 
On deeper quenches the system undergoes a continuous crossover to gel formation. We identify two regimes, 
'classical' phase separation, where single particle relaxation is faster than the dynamics of phase separation, 
and 'viscoelastic' phase separation, where demixing is slowed down appreciably due to slow dynamics in the 
colloid-rich phase. Particles at the surface of the strands of the network exhibit significantly greater mobility 
than those buried inside the gel strand which presents a method for coarsening. 



I. INTRODUCTION 

Dynamical arrest in soft matter can generally be divided 
into two classes, glasses and gels. Arrest is exhibited by a 
wide variety of systems, including those which can be treated 
as spheres with attractions which include colloid-polymer 
mixtures and protein solutions (2] [31. Dynamical arrest oc- 
curs upon either a global increase in density (vitrification) or 
arrested phase separation (gelation). In the latter case, phase 
separation is suppressed as the system undergoes arrest to 
form a non-equilibrium gel state, which we define as a per- 
colating assembly of particles (4|[5) which do not relax on the 
experimental timescale. 

In colloid-polymer mixtures, addition of sufficient polymer 
causes demixing into two or more phases (6l 0. In a one- 
component description, the degrees of freedom of the poly- 
mers are formally integrated out, leading to an effective at- 
traction between the colloidal particles whose range is fixed 
by the ratio of the polymer and colloid diameters #r and 
whose strength is set by the polymer concentration. To a rea- 
sonable approximation, polymer concentration plays the role 
of inverse temperature. The topology of the phase diagram is 
well-known: short-range interactions lead to only fluid-crystal 
phase coexistence, along with a metastable fluid-fluid binodal. 
However as the interaction range increases to #r w 0.3, fluid- 
fluid phase separation produces a stable colloidal liquid and 
colloidal gas CH3H1. 

While the topology of the phase diagram is well-known, 
the case of metastable amorphous states presents some chal- 
lenges to theory and simulation lUll . Smaller polymers lead to 
a shorter interaction range #r, so that upon demixing, the col- 
loids are closer together and thus the volume fraction of the 
colloid-rich phase c _ r is increased 021 [13]]. Deep quenches 



(higher polymer concentration) may also raise c _ r . During 
phase separation, if c _ r reaches volume fractions at which 
dynamical arrest occurs (0 c _ r ^ 0.57 — 0.59), demixing is sup- 
pressed and a gel results (UHUE). In the absence of crys- 
tallisation, such a system is nonetheless metastable to col- 
loidal liquid-gas coexistence. Fig. [T] shows a schematic sketch 
of the effect of interaction range on the phase behaviour (ig- 
noring crystal phases) l9lfT6l. 

Although such polydisperse dense (colloidal) liquids may 
be metastable to the formation of size-segregated crystals 
fTTlL to the best of our knowledge this has not been observed 
in experiments. Consequently, while most theoretical treat- 
ments necessarily consider the equilibrium monodisperse case 
where colloidal liquids are found up to c _ r < 0.5 (a regime 
in which dynamical arrest is not expected), the local vol- 
ume fraction in the 'arms' of a gel may be higher in experi- 
ments due to polydispersity. Although mean-field theories ex- 
ist where the metastable gas-liquid binodal and spinodal can 
be calculated lfl8lL such predictions have limited accuracy for 
many polymer-colloid size ratios iTTTl. 

Previous work has largely focused on gelation in systems 
with short-range attractions, typically with polymer-colloid 
size ratios #r < 0.2(4l [T9H231 such that there is no equilib- 
rium colloidal liquid phase. In such systems, gelation occurs 
upon crossing the metastable fluid-fluid spinodal 1 4 , l5l I2T1 l23l . 
In contrast, phase separation and its arrest in colloid-polymer 
mixtures with medium-to-long range attractions, where there 
is a stable colloidal liquid phase, are much less understood. 
Recent simulation (24) and experimental data(25l [26) suggest 
that a liquid-like colloid-rich phase exists, which grows con- 
tinuously in length scale as the system ages. 

For long-range (#r > 0.3) systems, we suggest two regimes 
of interest in relation to the gas-liquid binodal - shallow and 
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deep quenches. For shallow quenches, the kinetics of phase 
separation in colloid-polymer mixtures is analogous to liquid- 
liquid demixing in molecular systems l27ll . The structural 
evolution can be described by the Cahn-Hilliard theory of 
spinodal decomposition ESI E2l- Several regimes are dis- 
tinguished by the growth rate of the wavelength of density 
fluctuations. Immediately after quenching, a fastest-growing 
dominant mode q m appears with a wavelength Am, which ex- 
hibits time-dependent growth, with the power-law exponent 
changing from 1 /3 to 1 as the system evolves from the diffu- 
sive into the hydrodynamic coarsening regime. The charac- 
teristic domain size associated with the wavelength grows un- 
til it approaches the capillary length A c = (y/Apg) 1 / 2 , where 
7 is the interfacial tension and Ap the density difference be- 
tween the two phases, at which point any gravity-driven flow 
may rupture the bicontinuous structure (30). Previous work 
on spinodal decomposition in colloid-polymer mixtures has 
shown that the crossover between diffusive and hydrodynamic 
regimes is broad and continuous, and both theory and experi- 
ment indicate that, depending on the importance of hydrody- 
namic interactions, the exponent may vary between 0.2 and 

1.1 El eh. 

Quenching more deeply, the volume fraction of the col- 
loidal liquid increases |9l H8). The volume fraction of the 
colloid-rich phase may then approach (f) ~ 0.58, the point at 
which slow dynamics set in for hard spheres. Coarsening may 
therefore be suppressed and the bicontinuous network gener- 
ated by the initial spinodal decomposition becomes partially 
or, in the limit of very short-range systems, even completely 
arrested on the experimental timescale fT5l . The viscosity dif- 
ference between the phases may be reminiscent of the vis- 
coelastic phase separation proposed by Tanaka fT5i 
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FIG. 1 Schematic phase diagram showing the effect of attractive 
range on the gas-liquid coexistence line and the intersection between 
the coexistence line and the glassy region. In short-range systems 
(a), the liquid-gas binodal is quite flat and upon quenching, the 
colloid-rich phase rapidly approaches a regime of glassy dynamics 
1 12 1 . Conversely, in longer-range systems (b), the volume fraction 
of the colloid-rich phase increases more gradually upon quenching. 

In this work we explore these phase separation regimes in 
colloid-polymer mixtures with intermediate-to-long range at- 
tractions (0.45 <qR< 0.89) using confocal microscopy. Our 
data suggest that there is a continuous transition from gas- 
liquid phase separation to arrested phase separation (gelation) 
as a function of polymer concentration. The phase separation 
dynamics become increasingly slow as the system is quenched 
more deeply. Experimental phase diagrams show that as the 



attractive range increases, the region where phase separation 
is uninterrupted by gelation expands. In addition, the mech- 
anism of domain growth is studied in detail. We find that 
upon quenching, coarsening occurs predominately through 
enhanced particle mobility at the surface of the strands of the 
gel network. 



II. MATERIALS AND METHODS 
A. Colloid-polymer mixture I: Emulsions + HEC 

Our first model system consisted of a low polydispersity 
emulsion of poly(dimethyl siloxane) (Silicone DC200 fluid 
lOcSt, Serva) dispersed in a refractive index-matched mixed 
solvent of 1,2-ethane diol (EG) and water. The emulsion 
was stabilized against aggregation by a combination of three 
commercial non-ionic surfactants (w/w ratio 1:0.417:0.015); 
tristyrylphenol ethoxylate, 16EO (Tanatex), PEO-PPO-PEO, 
M w 4950, 30% PEO (Pluronic P103, BASF) and to give a 
small negative charge, the anionic surfactant sodium bis(2- 
ethyl 1-hexyl) sulfosuccinate (Aerosol OT, Cytec). Small 
droplets (PDMS: surfactant 1:0.12) had a hydrodynamic ra- 
dius g/2 = 142 nm and a size polydispersity of 0.10, while 
large droplets (PDMS: surfactant 1:0.05) had a hydrodynamic 
radius a/2 = 277 nm and a size polydispersity of 0.18. The 
index-matched continuous phase had a viscosity of 5.42 mPa s 
± 0.02 mPa s, as measured by capillary viscometry, a refrac- 
tive index of 1.396 at 543 nm, and a density of p c = 1.067 
g cm~ 3 ± 0.005 g cm -3 . The density of the PDMS oil is 
lower than the continuous phase (p m = 0.934 g cm -3 ± 0.006 
g cm -3 ) so that the dispersions cream in a gravitational field. 
The density mismatch of the oil is Ap = p c — p m = 0.13 ± 
0.01 g cm -3 . To screen out long-range electrostatic interac- 
tions and generate hard- sphere-like interactions between the 
emulsion drops we add 3 mM potassium chloride, which gives 
a Debye length of 6 nm. A dye, Sudan Black, was also 
added to increase contrast between colloid-rich and colloid- 
poor phases such that the phase separation could be identified 
by eye. Sudan Black preferentially adsorbs to the surface of 
the PDMS droplets, and so the colloid-rich phase appears dark 
and the colloid-poor phase light. 

Hydroxy ethylcellulose (HEC), which is a water-soluble 
non-ionic cellulose derivative, was used as the depletant. Two 
molecular weights, of hydroxy ethylcellulose were employed: 
Natrosol 250 M and Natrosol 250 HHX ( Ashland- Aqualon), 
with molecular weights of M w = 7.2 x 10 5 g mol -1 and 
M w = 1.3 x 10 6 g mol -1 respectively. Their estimated radii 
of gyration R g are 92 nm for Natrosol 250 M and 126 nm 
for Natrosol 250 HHX, as obtained by static light scattering. 
Small droplets combined with Natrosol HHX gives a polymer- 
colloid size ratio, #r = (2R g /a), of 0.89; small droplets with 
Natrosol M gives #r = 0.65; finally, large droplets with Na- 
trosol HHX gives q R = 0.45. 

Samples were prepared by first diluting a stock solution of 
the polymer with solvent, ensuring proper dispersion of the 
polymer. Then, solvent containing fluorescent rhodamine-B 
dye (Sigma Aldrich) was added such that the concentration of 
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rhodamine-B was 0.02 gL _1 in the aqueous phase for all sam- 
ples. Finally, stock emulsion was added and the sample gently 
stirred for at least ten minutes before being quickly transferred 
to the imaging cell. The start of the experiment (t = 0) was 
taken as the instant when stirring ceased. 



B. Colloid-polymer mixture II: PMMA + PS 

For our second system, we used poly(methyl methacrylate) 
(PMMA) core-shell particles (32||33, sterically stabilised with 
poly(12-hydroxy stearic acid) (PHSA), with a mean radius 
416 nm and a polydispersity of 0.10, as determined by scan- 
ning electron microscopy. Both the core and shell were la- 
beled with fluorescent dyes: the core with coumarin and the 
shell with rhodamine-B. The particles were suspended in a 
density-matched 78/22 (%w/w) mixture of cyclohexyl bro- 
mide and cis-decalin. 4 mM of tetrabutyl ammonium bro- 
mide was dissolved in the solvent mixture to screen electro- 
static interactions and ensure hard-sphere-like behaviour of 
the colloids. Depletion attractions were induced by addition 
of polystyrene (PS) as a non-adsorbing polymer. The polymer 
molecular weight was 2.06 x 10 7 g mol -1 . We estimate the 
polymer radius of gyration, R g , to be 203 nm[34], such that 
the polymer-colloid size ratio #r was 0.49, which was compa- 
rable to that of the smallest size ratio in the emulsion droplet 
system. 



C. Confocal microscopy 

Confocal microscopy was used to study the temporal evolu- 
tion of the bicontinuous network formed in our aqueous emul- 
sion samples. A supporting frame was constructed to allow a 
light microscope (Zeiss, Axioskop SI 00) fitted with a confo- 
cal scanning head (Zeiss, LSM Pascal) to be mounted hori- 
zontally. A laser wavelength of 532 nm was used to excite 
rhodamine molecules dissolved in the solvent, such that the 
background appears light and the emulsion droplets dark. For 
clarity, we have inverted the micrographs pertaining to system 
I so that the colloid-rich phase appears bright and the colloid- 
poor phase dark. Emulsion-polymer dispersions (system I) 
were studied in large rectangular glass cuvettes with the front 
wall of the cell constructed from a 170 jim thick optical qual- 
ity cover glass. The sealed cuvette had a square cross- section 
with an internal dimension of 13 mm, a height of 30 mm, and 
contained a small hole at the top of the cell through which the 
cell was filled. 

For PMMA samples (system II), a conventionally mounted 
confocal microscope (Leica SP5 with resonant scanner) was 
used to image in the xy-plane and a laser with wavelength 
488 nm was used to excite fluorescence from the coumarin- 
labelled core of the particles. 2D Images were collected at 
time intervals of St =150-300 milliseconds over a total time 
period of 10-20 seconds, with typically 100-400 particles in 
each frame. Data was averaged over 10-20 runs in order to 
improve statistics. 



D. Analysis 

To study the structural evolution quantitatively, length 
scales characteristic of the gel structure were extracted from 
the 2D-confocal images. For an image of width L taken at time 
t, we calculated the radially-averaged structure factor S(q,t) 

%,') = / dq , (/(q , ,0/(-q , ,0> (D 

2nqAq J q <\ q >\< q +A q 

where /(q,f) is the 2D Fourier Transform (FT) of the image 
intensity /(r,f) and Aq = 2n/L. The 2D Fourier transform 
displays a ring of high intensity at low q which signifies the 
presence of a large characteristic length scale. We characterise 
the dominant wavevector by the first moment (q(t)) of the 
structure factor, where 

(q(t)) = [ b dqqS(q,t)/ t &qS{q,t) (2) 

J a J a 

with limits a = 0.02 /im -1 and b = 0.4 jUm -1 . The domi- 
nant wavevector was then converted into a characteristic size 
X c (t)=2K/{q(t)). 

In the PMMA system, the motion of the colloidal parti- 
cles was tracked using modified versions of standard tracking 
routines l35l l36l . Particles were separated into two groups: 
we define surface particles as those with n = 1 — 4 neigh- 
bours throughout the image set, and bulk particles as ones with 
n > 4. 

In order for timescales to be compared directly in sys- 
tems with different particle sizes and solvent viscosities, we 
quote all times in units of the characteristic Brownian time 
Tb = 37T77LCT 3 / (4k#T), which is comparable to the time taken 
for a particle to diffuse its own radius in a solvent of viscos- 
ity 7]l. For the PDMS emulsion system, the long time vis- 
cosity 7]l of the continuous phase was determined from the 
diffusion constant of particles of known radius suspended in 
HEC solutions, using dynamic light scattering. Viscosities T]l 
were deduced from the measured diffusion constant using the 
Stokes-Einstein equation D = k^T / (37irf^G) . The concentra- 
tion dependence of the viscosity was then fitted to the Martin 
equation 03, tj l /t]o = 1 + [ri]C v exp(kn[ri}C p ), where tj l is 
the limiting low shear viscosity, T]o is the solvent viscosity, 
[rj] is the intrinsic viscosity and ku is the Huggins constant. 
This equation is a reasonable prediction of polymer viscos- 
ity in the semi-dilute regime. Regression gave fitted parame- 
ters of [tj] = 0.30 Lg _1 and kn = 1.2 for the HHX grade, and 
[tj] = 0.25 Lg _1 , and ku = 1 .2 for the lower molecular weight 
M grade. The estimated Brownian times Tb were between 0.1 
s and 1.1s, depending on the particle size, polymer molecular 
weight and polymer concentration. In the PMMA system, the 
Brownian time was determined as Tb =0.7 s. 

E. Phase diagrams 

Experimental phase diagrams at size ratios #r = 0.45 and 
#r = 0.89 were constructed for system I. The overall col- 
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loid volume fraction in each sample was fixed at = 0.2 for 
all samples, and polymer concentration C p and hence quench 
depth was varied, such that all state points lie on a vertical line 
through the phase diagram. Phase diagrams were obtained by 
allowing emulsion droplet samples shallowly quenched into 
the two-phase region to fully equilibrate, then measuring the 
colloid volume fraction in both phases using UV- visible spec- 
troscopy to measure absorbance of Sudan Black at 600 nm, 
which was proportional to the droplet concentration. Fig. [2] 
shows the phase diagrams we obtained. Note that we fix the 
overall (f) in each sample. Thus our representation is approxi- 
mately analogous to the polymer reservoir representation (9j, 
scaled by the free volume fraction for (j) =0.2. 

The gas-liquid coexistence boundaries at = 0.2 were at 
C p = 0.40 gh~ l and C p = 0.60 gL~\ for q R = 0.45 and q R = 
0.89 respectively. As the time taken to equilibrate increases 
rapidly with polymer concentration, only shallow quenches 
were experimentally accessible for this method and we note 
that underestimation of the colloidal volume fraction of the 
colloid-rich phase cannot be ruled out. Although the polymer 
concentrations used here are less than those used for the con- 
focal experiments, extrapolations are sketched to show that 
the in the case of the smaller q, colloid-rich phase becomes 
dense with relatively shallow quenches after the coexistence 
line is crossed. In contrast, the dense branch of the long-range 
system requires much deeper quenches to achieve comparable 
volume fraction. 

Confocal microscopy was used to identify state points as ei- 
ther one- or two-phase and therefore to determine an estimate 
of the phase boundary for q R = 0.65 (system I) and also for 
system II. The phase boundary for q R = 0.65 was estimated 
as C p = 1.00 ±0.20 gL _1 at (j) = 0.2. Now the phase bound- 
ary for q R = 0.65 in units of gL -1 does not lie between those 
corresponding q R = 0.45 and q R = 0.89. This is due to the 
different polymer molecular weights used to obtain these size 
ratios. Similarly for system II, the phase boundary was esti- 
mated as C p = 0.13 ±0.03 gL" 1 . 

As the molecular weights and the solvency conditions of 
the polymers varied depending on q R and also on the system 
used, polymer concentrations are scaled by the polymer con- 
centration at the (colloidal) gas-liquid critical point C cr i t , so 
as to allow direct comparisons to be made. We estimate C cr i t 
by taking the polymer concentration at the liquid-gas phase 
boundary at (j) = 0.2. Since our phase diagrams (Fig. [2]) are 
plotted analogously to the polymer reservoir represenation, 
close to criticality the phase boundary is quite flat with re- 
spect to (j) so that (j) = 0.2 provides a reasonable estimate for 
the critical volume fraction cr i t [ 38, 39 ]. For system I, C cr i t = 
0.40 gL" 1 , 1.00 gL- 1 and 0.60 gL" 1 for q R = 0.49, q R = 0.6 
and q R = 0.89 respectively. For system II, C cr i t = 0.13 gL -1 . 



III. RESULTS AND DISCUSSION 



of system I (not density-matched), instead of gel formation, 
gravity-driven flow occurs almost immediately after a bicon- 
tinuous network forms, reminiscent of images showing the 
Rayleigh-Taylor instability in colloidal suspensions l30l l40l . 
Deeper quenches allow a stress-supporting network to form, 
where gel-like structures evolved ~30 seconds after the sam- 
ple mixing is ceased. Shorter-range, strongly attractive sam- 
ples produced a long-lived network of fine strands, whereas 
samples with longer-range attractions showed a continuously 
changing viscoelastic gel structure. 

At a polymer concentration of C p /C cr j t = 1.20 (Fig. [5]) 
in the density-matched system II, liquid-gas phase separation 
progresses to completion over about two hours, with a per- 
colating network of the colloid-rich phase coexisting with a 
colloidal gas, highly reminiscent of spinodal-like patterns also 
seen in the emulsion droplet system (I). 

The rate of structural coarsening in system I can be rep- 
resented by the growth of the characteristic domain size, as 
defined by the dominant wavevector of the 2D Fourier Trans- 
form of each confocal image. Each data set exhibited a pe- 
riod of power law growth, where the rate of growth was given 
by the power law exponent. Fig. [4] shows examples of how 
growth exponents were determined. 

Data deviated from power law growth at very short times 
(< 100tb), where the domain size remained nearly constant. 
This is indicative of a linear Cahn regime in early- stage spin- 
odal decomposition and has previously been observed in other 
colloidal systems l26ll4T1l . 

The time-dependence of domain size A c for varying poly- 
mer concentrations and size ratios is shown in Fig. [5] Growth 
exponents as a function of polymer concentration are shown 
in Fig. [6] At low polymer concentration, a regime of fluid- 
fluid phase separation was observed for all size ratios. In 
this regime, the growth exponent is independent of quench 
depth for a given size ratio but ranged between 0.2 and 0.6 
for the range of q R studied. This is consistent with previ- 
ous work on phase separation in colloid-polymer mixtures 
where the growth exponent varied between 0.2 and 1.1 due 
to the broad crossover between diffusive and hydrodynamic 
regimes l27l[3T1l . 

B. Time evolution: connection to phase behaviour 

To gain insight into the phase separation dynamics, we as- 
sume these are related to the structural relaxation time in the 
colloid-rich phase. Neglecting the interfaces resulting from 
the bicontinuous network, the colloid-rich phase may be ap- 
proximated by hard spheres at the same volume fraction, as 
the polymer concentration in this phase is low 1 9 ] . The struc- 
tural relaxation time T a in the colloid-rich phase can then 
be estimated, as a function of the volume fraction 0, by the 
Vogel-Fulcher-Tammann equation 1 42T- I441 : 



A. Time evolution: observation 

At short times (>1 minute), spinodal decomposition pat- 
terns were observed. For very shallow quenches, in the case 
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Equation [3] is plotted in the inset in Fig. fusing parameters 



FIG. 2 Phase diagrams at constant overall sample colloidal volume fraction (j) = 0.2, (a)#R = 0.45 and (b)#R = 0.89. In each phase diagram, 
the left-hand y-axis gives polymer concentrations rescaled in terms of the critical polymer concentration C cr i t as estimated from the phase 
diagram, and the right-hand y-axis gives polymer concentration in units of gL _1 . Dotted lines show qualitative extrapolations of the phase 
boundary for deeper quenches. Filled-in data points indicate data from phase separation experiments, unfilled data points show state points 
used for confocal imaging and Fourier lengthscale data. Inset shows dependence of structural relaxation time z a on colloid volume density 
for hard spheres. Shaded region in inset denotes the range of = 0.57 — 0.59 where the dynamics become very slow. 
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FIG. 3 PMMA particles undergoing fluid-fluid phase separation, 
Cp/Qrit = 1-20, = 0.2, #r = 0.49 showing characteristic spinodal 
pattern, with inset showing detail. Scale bar in top left-hand corner 
indicates 25 ftm; scale bar in inset indicates 5 jim. 



from ref. El, namely A = 2.25 x 1(T 2 , 8 2.0 and O 
0.637. As depicted in the inset to Fig. |2jb), when (j) ~ 0.57 — 
0.59, relaxation times increase very strongly. 

The dashed lines in Figs. [2] (a) and (b) indicate the expected 
increase in volume fraction of the colloid-rich phase c _ r at the 
state points considered in our time-evolution (Fig. [5]). Our ex- 
trapolation suggests that, in the intermediate and slow phase 
separation regimes, the volume fraction of the colloid-rich 
phase reaches c _ r ~ 0.57 — 0.59, such that the dynamics of 
phase separation is limited by particle relaxation in the dense 
phase. In other words, the particle level dynamics become 



FIG. 4 Examples of growth exponent fits for selected domain size 
growth datasets. Here we plot the characteristic domain size Xq as a 
function of time. #r = 0.65, (j) = 0.2. Growth exponents of 
0.229±0.011, 0.132±0.006, 0.060±0.004 respectively. 



coupled to phase separation as a whole. Conversely, at weaker 
quenches where the volume fraction of the colloid-rich phase 
is not high enough to exhibit slow dynamics, T a is very much 
less than the dynamics of phase separation. In short, our sys- 
tem exhibits two regimes of phase separation - 'normal' and 
'viscoelastic' l45ll . 

This crossover between phase separation and gelation gives 
rise to an intermediate regime where phase separation still 
progresses to completion, but in contrast to ordinary fluid- 
fluid phase separation, the growth exponent decreases with 
increasing polymer concentration. This regime is dominated 
by slow dynamics in the colloid-rich phase, where increasing 
densities result in a viscous yet still fluid state. A correspond- 
ing increase in T a can be seen, although in this regime T a is 
still easily accessible on the experimental timescale. These 
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regimes are indicated in Fig. [6] 

The onset of gelation can be identified from the data by 
a decrease in the growth exponent during the imaging time, 
at t ~ 100 seconds (which corresponds to around IOOOtb) 
as shown in Fig. [5] In these samples, phase separation 
is suppressed by gelation after the colloid-rich phase densi- 
ties to form a glass. This type of behaviour has also been 
seen in simulations! 24 ] of the Lennard- Jones system, as well 
as experiments on systems with short-range attractions l4l. 
Quenching more deeply, the slow growth regime of character- 
istic domain size proceeds is quickly reaching after quenching 
and so the change in growth exponent becomes too fast to be 
experimentally accessible. 

In the gel regime, the structure continues to coarsen, i.e. 
the growth exponent is non-zero. This differs markedly from 
the behaviour of gels with short-range attractions and, at first 
glance, may appear paradoxical, as the gel structure is com- 
posed of a glassy network where particles are unable to rear- 
range. 



C. Local dynamics 

We now address this apparent paradox by studying the mo- 
tion of individual particles using experimental system II. We 
find that that particle mobility at the surface of the gel strand 
is enhanced in comparison to particles in the interior of the gel 
strand, such that occasional particle-hopping was observed at 
the surface of the gel strands even in gels that were 'arrested'. 

Dynamical heterogeneities became significant as the sys- 
tem approached gelation. At a polymer concentration of 
Cp/Qrit = 2.40, it was apparent that the system consisted of 
two particle populations, an effect previously observed in col- 
loidal systems l23l l46l l47l . and also for polymer glasses [48]. 
'Slow' particles were stuck in an arrested network, and 'fast' 
particles, which were able to diffuse over a number of parti- 
cle diameters, occupied surface sites on this network. Mean 
square displacements of bulk and surface particles (Fig. [7]) 
show that they correspond to 'slow' and 'fast' particles re- 
spectively. Sub-diffusive behaviour is observed for all parti- 
cles, although this may be due to the relatively short imaging 
time. Due to the experimental timescale, the extent of parti- 
cle diffusion is limited, but there are indications that a plateau 
may be developing in the data for C p /C cr i t = 2.40. As quench 
depth increased, the number of 'fast' particles decreased. At a 
deep quench of C p /C cr i t = 9.70, no fast particles were detected 
over a two-hour period. 

Some localised areas of the gel network surface exhibited 
greater particle mobility than others, and particles appear to 
rearrange cooperatively, an example of which is shown in Fig. 
[8] It is possible that self-generated mechanical stress within 
the gel and the resulting deformation lowers the activation en- 
ergy of particle-hopping in localised areas. It has previously 
been proposed [ 49 ] that self-generated stress builds up at local 
weak points within a transient gel network, thereby increasing 
particle mobility in small regions and providing a mechanism 
for gel collapse in arrested gels. 
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FIG. 5 Here we plot the characteristic domain size A c , obtained 
from 2D-FFT of confocal images of system I, as a function of time. 
= 0.2, (a) q R = 0.45, (b) q R = 0.65, (c) q R = 0.89. Indicated in 
(a-c) are estimated volume fractions in the colloid-rich phase c _ r , as 
compared volume fraction at which hard spheres undergo the glass 
transition, g . 



IV. CONCLUSIONS 

We have presented an experimental study of dynamics in 
phase- separating colloid-polymer mixtures with interaction 
range between 0.45 and 0.89. Our data suggest that the 
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FIG. 6 Growth exponent of domain sizes as a function of polymer 
concentration, #r = 0.45 — 0.89, (j) = 0.2. Dashed lines are a guide 
to the eye. Labels indicate behavioural regimes as follows: (a) 
Normal phase separation, (b)Intermediate 'slow' phase separation 
and (c) gelation. 
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FIG. 8 Localised particle-hopping on surface of arrested gel 
network. Magnified snapshots of an area (small box) are shown in 
the four larger boxes as a time-sequence. The particles undergo a 
local re- arrangement over ~ 500 s, as indicated by the outlined 
particles. The snapshots are taken at t = 320 s, 410 s, 420 s and 460 
s after imaging began, which correspond to 460tb, 590tb, 600tb 
and 660t b . C p /C crit = 3.90, q R = 0.49, = 0.2. Scale bar in top 
left-hand corner indicates 10 jUm. 



FIG. 7 Mean square displacements for C p /C cr it = 1.20 and 
Cp/Qrit = 2-40 in units of jUm 2 and particle diameter a, with 
particles grouped according to number of neighbours n. Dotted line 
shows diffusive linear growth of MSD. = 0.2, q R = 0.49 



crossover from phase separation to gelation is a continuous 
process as function of polymer concentration, with phase sep- 
aration dynamics slowing gradually until arrest is reached. By 
analysing the time growth of the domain size during phase 
separation, we observe a range of behaviour from liquid-gas 
phase separation at shallow quench depths to dynamic arrest 
and gelation at deep quenches. Data from confocal imaging is 
supported by experimental equilibrium phase diagrams, in or- 
der to highlight the changes in dynamics with attractive range. 
As the attractive range increases, the region of the phase dia- 
gram corresponding to phase separation becomes larger, and 
the crossover into gelation more gradual. 

In addition, we find evidence of a mechanism of structural 
coarsening and domain growth in 'arrested' gels, where re- 
arrangements on the surface of the gel strand may continue 
despite a glassy, arrested state inside the gel strand, an ef- 
fect which becomes more significant with increasing attractive 
range. 

The dynamics of the colloid-rich phase are taken to be simi- 



lar to hard spheres at the same estimated volume fraction. Un- 
der this assumption we find the dynamics of phase separation 
couple to the structural relaxation at the particle level in the 
colloid-rich phase upon deep quenches where slow dynamics 
set in. At weaker quenches, the phase separation dynamics 
decouple from the single-particle dynamics, and in this way 
we distinguish two modes of phase separation - 'normal' and 
'viscoelastic'. Further work could investigate this observation 
more closely. In particular, a more accurate determination of 
the volume fraction of the colloid-rich phase would be desir- 
able. While this may be tackled experimentally or by simula- 
tion, theoretical extensions of metastable liquid binodals and 
spinodals would be most helpful. While mean-field theory 
fT8l is readily available for the range of volume fractions in 
the colloid-rich phase relevant to gelation (0 c _ r > 0.57 — 0.59), 
especially for smaller size ratios, its accuracy may be limited 
full . We hope to have raised the issue that colloid-polymer 
mixtures (and other model systems such as DNA-coated col- 
loids [50] frequently become involved in metastable liquid 
states at high volume fraction, and that treatment of such states 
with theories more sophisticated than mean-field (such as hi- 
erarchical reference theory |51|) or self-consistent Ornstein- 
Zernike approximation [ 13 ] may prove invaluable if extended 
to higher colloid volume fraction for a variety of interction 
ranges. 
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